A bosonization approach for bilayer quantum Hall systems at vt = 1 
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We develop a non-perturbative bosonization approach for bilayer quantum Hall systems at z/t = 1, 
which allows us to systematically study the existence of an exciton condensate in these systems. An 
effective boson model is derived and the excitation spectrum is calculated both in the Bogoliubov 
and in the Popov approximations. In the latter case, we show that the ground state of the system 
is an exciton condensate only when the distance between the layers is very small compared to the 
magnetic length, indicating that the system possibly undergoes another phase transition before the 
incompressible-compressible one. The effect of a finite electron interlayer tunnelling is included and 
a quantitative phase diagram is proposed. 
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Introduction. The recent conjecture of Bose-Einstein 
condensation (BEC) of excitons in a bilayer quantum 
Hall system (QHS) at total filling factor one {vt = 1) 
has attracted a great deal of attention U Q • This sys- 
tem consists of two quantum wells (layers) separated by 
a fixed distance d, under a perpendicular magnetic field 
B. Each layer has filling factor v = n(j>o/B = 1/2, where 
n is the layer electron density and (/>o ~ he/ e is the mag- 
netic flux quantum. By increasing the ratio d/l, where 
I = yJhc/eB is the magnetic length, the bilayer QHS at 
Vt — \ undergoes a transition from an incompressible to 
a compressible phase 0. Indeed, for small d/l, sponta- 
neous interlayer coherence develops and the system be- 
haves as a single-layer QHS with v — \, whereas for large 
values of d/l intralayer correlations become important 
and the system behaves as two independent QHSs with 

The nature of the ground state of the bilayer QHS at 
i^T = 1 is an issue not completely settled yet. From the 
experimental point of view, tunnelling conductance ex- 
periments 1^ show a zero bias peak which does not seem 
to diverge with decreasing temperatures. In addition, a 
linear current-voltage characteristic, instead of a power 
law one, has been observed in magneto- transport mea- 
surements 0, 0] ■ Both points raise doubts about the ex- 
istence of a true exciton condensate, which might have su- 
perfluid properties 0. Theoretically, diagrammatic cal- 
culations which rely on the existence of an exciton 
condensate, indicate that the ground state of the system 
becomes unstable for d/l > 1.2. Single- mode approxima- 
tion 10] and generalized random phase approximation 
calculations JJJ confirm the result. Recently, Fertig and 
Murthy suggested that this "imperfect" two-dimensional 
superfluid behavior of the bilayer may be understood in 
terms of a coherence network induced by disorder, which 
breaks up the system into large and small regions with, 
respectively, weak and strong interlayer coherence [l^. 

In this letter, we develop a non-perturbative bosoniza- 



tion formalism for bilayer QHSs at vt — ^, which allows 
us to properly explore the idea of exciton condensation in 
this system. An interacting boson model is derived and 
the excitation spectrum is calculated both in the Bogoli- 
ubov and in the Popov approximations. In the latter 
case, we find that the ground state of the system is a 
true exciton condensate only when d/l < {d/l)c. For 
the zero electron interlayer tunnelling case, {d/l)c — 0.4, 
but this parameter increases with the tunnelling strength. 
This new phase transition (or crossover) takes place at 
a critical ratio {d/l)c much smaller than {d/l)]r'~^, where 
the incompressible-compressible phase transition is ex- 
perimentally observed. Our findings, which are based on 
a proper treatment of the Coulomb interaction, yield a 
quantitative phase diagram for the bilayers. 

The model. Let us consider a bilayer system with N 
spinless electrons moving in the {x,y,z = 0) plane and 
N spinless electrons moving in the {x, y, z — d) plane 
in an external magnetic field B = Bz. For each layer, 
we restrict the Hilbert space to the lowest Landau level 
(LLL) and consider TV = N0/2, where is the Landau 
level degeneracy. We also introduce a pseudospin index 
(tj =t, i) in order to define to which layer each electron 
is associated. 

The Hamiltonian of the system is given by 



H = Ht + V, 



(1) 



where Tix describes the electron tunnelling between the 
two layers. 



(2) 
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and V is the Coulomb interaction term (unit area system) 
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Here Ag^s is the electron interlayer tunnelling, the 
fermion operator cj,j ^ creates an electron with guiding 
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center m in the LLL of the a layer (see Fig.^, Pa{k) is 
the Fourier transform of the projected density operator of 
pseudospin a electrons, and Wff(fc) = '^u(^) ~ VA{k) — 
(27re2/eA:)exp(-|fcZ|72) and v^i{k) = v^{k) = VE{k) = 
(27re^/e/c) exp(— |fcZp/2) exp(— fed) are, respectively, the 
Fourier transform of the intrAlayer and intErlayer inter- 
action potential. 

Firstly, we will concentrate on the limit Ag^^s = 0. 
The Hamiltonian reduces to the Coulomb term Q 
which can be rewritten in terms of the total electron 
density p{k) = p-\{k) + pi{k) and the z-component of the 
pseudospin density Sz{k) — [p]{k) — Pi{k)\ /2 operators 
as 

^ = \Y. Mk)p{k)p{-k) + 2 J2 vz{k)Sz{k)Sz{-k), 

k k 

(4) 

with 

Vo/z{k) = \[vA{k) ± VE{k)] = !^e-l'='l'/2(l ± e-'^"). 

The restriction of the Hilbert space to the LLL together 
with the introduction of the pseudospin language renders 
the description of the bilayer QHS at i^t = 1 analogous 
to the one of the single-layer QHS aXv = \ when the elec- 
tron spin degree of freedom is included. For the latter, it 
was shown that the particle-hole pair excitations (mag- 
netic excitons) of the ground state (quantum Hall ferro- 
magnet) can be approximately treated as bosons jl."^ . In 
this framework, the electron and the z-component of the 
(pseudo)spin density operators are written as 

p{k) = 4,oN0 + 2i^sin(kAq/2)6j,+q6q, 
q 

Sz{k) = i4,oN^-5Ic°^('^^q/2)^k+q^' (5) 
q 

with k A q = Pz • (k x q). Here, 6^ and 6q are boson 
operators which obey the canonical commutation rela- 
tions [6^,6q,] = [6q,6q'] = 0, and [6q,5^,] = 5c^,c^'■ When 
6^ is applied to the quantum Hall ferromagnet, it cre- 
ates a magnetic exciton whose momentum q is related to 
the vector r between the guiding centers of the excited 
electron and hole by |r| = f'\ci\ [Fig. [ija)] In this 

formalism, the configuration of a bilayer QHS aX vt = ^ 
corresponds to a system with N0/2 bosons. 

Substituting Eqs. (O into Q and normal-ordering the 
result, we obtain, apart from a constant related to the 
positive background, an interacting boson model 

Hb = ^Wqh\b^ + ^ 2wk(p,q)&],+p5|j_i^6q6p,(6) 

q k,p,q 

where the dispersion relation of the bosons is given by 
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FIG. 1; Schematic representation of (a) one-boson state (mag- 
netic exciton) of the quantum Hall ferromagnet with |qZ^| — 1 
and (b) the condensate of zero-momentum bosons. The 

guiding center quantum number is denoted by m. 

with Jo(a;) denoting the Bessel function of first kind. The 
second term of Wq is similar to the capacitor-like one 
introduced phenomenologically by Girvin ^] . The boson 
interaction potential is given by 

^'k(p, q) = wo(fc) sin(k A p/2) sin(k A q/2) 

+ wz(fc)cos(kAp/2)cos(kAq/2). (7) 

By taking the limit d/l ^ in Eq. ((HJ, we recover the 
boson model for the single-layer QHS at v = 1 01 • 

Bogoliubov approximation. We start the analysis of the 
interacting boson model © within the Bogoliubov ap- 
proximation We assume that the bosons condense 
in their lowest energy state, the q = mode, which means 
that the ground state of the system is roughly given by 
the state shown in Fig.^b). This assumption is in agree- 
ment with the scenario proposed in Refs. 2J. Therefore, 
the boson operators fiq^g and 6q=o can be replaced by 

complex numbers, i.e., 6q,6o — » VNq, where Nq is the 
(macroscopic) number of bosons in the q = mode. 

After substituting 6q and bo by v^Nq in the Hamil- 
tonian jni, including the chemical potential p explicitly, 
i.e., Hb — > K = TiB—pN, and keeping only the quadratic 
terms in the boson operators, we obtain 

K « 2t;z(0)N2 + (u>o-Ai)No-(l/2)^e, 

+ 2 ^ [^9(^q^q + ^-q^-q) + -^«(^l^-q + ^-q^q)] - 

(8) 

where Cq = Wq — p + 4Not'z(0) + 4NoWz(g) and A, — 
4NoWz(9). The above Hamiltonian can be diagonal- 
ized by the canonical transformation 6q — cosh(7q)aq — 
sinh(7q)a^lq, and 6^ = cosh(7g)aJj — sinh(7q)a_q, with 7^ 
real, and therefore the dispersion relation of the quasi- 
particles is given by 
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FIG. 2; (color online) Dispersion relation of the quasiparticles 
in the Bogoliubov approximation for different values of the 
ratio d/l: 0.2 (black), 0.5 (red), 0.8 (green), 1.0 (blue), 1.5 
(violet), and 2.0 (magenta), (a) Asas = and (b) Asas = 
0.1 e^/d. 

In the one- loop approximation (l^, /i = wq + 4v2(0)No, 
and thus, eq — Wq — wq + 4NoUz('?). Using that {N) = 
No + X]q/o(^q^q) ~ N0/2, we can calculate Nq, and de- 
termine rig. The latter is illustrated in Fig.|2Ia) for sev- 
eral values of d/l. 

In the long wavelength limit, Q.q ~ hvq, where the 
linear mode velocity is v — y^2noeBd/l{l/h)e'^ /el. Here 
es = -d/4l+y^^i/32{l+d^/P)eicp {d^/2P) Erfc(d/\/2Z), 
and Erfc(a;) is the complementary error function. For 
d/l = 0.2, 0.5, 0.8, 1.0, 1.5, and 2.0, we find v = 4.47, 
4.64, 3.92, 3.39, 2.27, and 1.52 x 10"* m/s, respectively. 

The dispersion relation Qq of the quasiparticles in- 
creases with momentum for all values of d//, in qualita- 
tive agreement with calculations done in Refs. |9l ll0llll| 
in the region d/l < 1.0. However, these previous cal- 
culations found that a minimum (roton-like excitation) 
appears around |q/| ^ 1.0 for d/l > 1.0 and that this 
minimum becomes soft at d/l > 1.2, whereas, in our 
case, no minimum arises when d/l increases. The behav- 
ior of the condensate fraction uq — No/(N0/2) in terms 
of d/l shows that the Bogoliubov approximation for the 
boson model JHJ holds only for small d/l. As the ratio 
d/l increases, the number of bosons in the condensate 
decreases continuously (no < 0.3 for d/l > 1.0), indi- 
cating that the interaction between excited (out of the 
condensate) bosons starts to play an important role in 



the description of the system. 

Popov approximation. It is possible to go a step further 
in the analysis of the interacting boson Hamiltonian |(BJ , 
by including the interaction between the excited bosons, 
which is neglected in the Bogoliubov approximation. 

Let us consider the four-operator terms of the Hamil- 
tonian © with |q^| ^ and treat them on the mean- 
field level, including only normal averages. This treat- 
ment resembles the so-called first-order Popov approx- 
imation for Bose gases at finite temperatures [l2J. By 
applying this approximation to the boson model ®, we 
obtain a Hamiltonian similar to Eq. JHJ with the re- 
placement €q ^ Iq = Wq - ^ + A'HqVz{0) + A^QVz{q) + 

4Ep5^q (^^^(0) + Vp-c^{^.■, q)) {hybp). After diagonalizing 
the new Hamiltonian and calculating the chemical poten- 
tial, we find the following set of self-consistent equations 

£9 = Wq- Wo + 2'N^vz{q) 

+ 4^bp-q(q,q)-wz(p)-wz(g)]^p, 

= 7^1 -A^ (10) 

Tq = {eq/^q-\)/2, 

where Clq is the (new) quasiparticle dispersion relation. 

Solutions for the above self-consistent problem can only 
be obtained for d/l < {d/l)c — 0.4. In this case, the 
dispersion relation of the quasiparticles also vanishes lin- 
early with momentum for all d/l, qualitatively confirming 
the results obtained within the Bogoliubov approxima- 
tion. The velocities of the linear modes are renormalized, 
for instance, v = 4.30 x 10^ m/s for d/l = 0.2. 

However, as the ratio d/l increases, it is not possi- 
ble to find real solutions for Eqs. 1)10(1 . implying that 
the ground state of the system is no longer a true con- 
densate of bosons with zero momentum. This indi- 
cates that a new phase sets in above {d/l)c and below 
the well-known incompressible-compressible phase tran- 
sition, which experimentally takes place at {d/l)l~'-' . No- 
tice that {d/l)c = 0.4 is much smaller than the critical 
ratio {d/l)l~'~' , which, in samples with negligible electron 
interlayer tunnelling, is around 1.6 — 1.8 !al3|i indicat- 
ing that the novel phase should be observable within a 
large region in the parameter space. The ground state of 
the bilayer QHS dX vt — ^ in the region {d/l)c < d/l < 
{d/l)l~'^ should then be more complex than a pure ex- 
citon condensate, possibly a zero momentum boson con- 
densate coexistin g w ith a charge-density wave state as 
suggested in Ref. [l8| . A detailed study of the properties 
of this novel phase based on the interacting boson model 
© and the analysis of the nature of this new phase tran- 
sition (or crossover) will be the subject of further inves- 
tigations. 
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FIG. 3: Phase diagram (djl x Asas) of the bilayer QHS 
at vt = the filled circles are the calcultated critical ratio 
{d/l)c (the soUd line is a guide to the eyes) and the dashed 
line is the experimental estimate for {d/l)i~^ (from Ref. Q)- 



Our findings open up the possibility of the existence 
of a two-component phase in bilayers, based solely in an 
appro pria te treatment of the Coulomb interaction. Dis- 
order though probably relevant, does not necessarily 
need to be invoked in order to generate a richer phase di- 
agram. 

Finally, we should mention that our findings are quite 
different from previous theoretical calculations [sl, [ifll, lll| 
which point out that an exciton condensate is stable for 
d/l < 1.2. Within our approach, this instability app ears 
at even smaller values of d/l. Moreover, in Ref. [10|, the 
softness of the roton-like excitation at d/l ~ 1.2 is asso- 
ciated with the incompressible-compressible phase tran- 
sition, which is not the case here. 

Finite electron interlayer tunnelling. The effect of a fi- 
nite electron interlayer tunnelling (Ag^g ^ 0) can be eas- 
ily included in our formalism. Now the complete Hamil- 
tonian should be considered. It is easy to show that 
the bosonic representation of Ht [Eq. ©] is quadratic 
in the boson operators. Adding this extra term to the 
boson model and following the same steps as before, 
we find, in the Bogoliubov approximation, that the dis- 
persion relation of the quasiparticles is also given by Eq. 
®, but now €q = A5^s(l/4-hl/no)vW72 + Wq-wo + 
4J>ioVz{q). Notice that the spectrum is no longer gapless. 
As AsAS increases, a minimum (roton-like excitation) 
develops in the excitation spectrum around |q^| = 2 for 
d/l > 1.0. The energy of the minimum reduces as the 
ratio d/l increases. The case Ag^s = 0.1 e'^/el is illus- 
trated in Fig.|2Ib) for several values of d/l. 

In the Popov approximation, the final set of self- 
consistent equations are similar to Eqs. 1)10(1 with the 
replacement Cq ^ Iq ~ iq + As^s(3no -I- 2)/\/32no. The 
inclusion of a finite eletron interlayer tunnelling increases 
the critical ratio {d/l)c (see Fig. OJ. However, the novel 
phase should remain robust in a sizable region because 
{d/l)l~'~^ also increases with the tunnelling strength as 
verified experimentally 3]. For d/l < {d/l)c, the dis- 



persion relation of the quasiparticles is quite similar to 
the results found in the Bogoliubov approximation apart 
from small renormalizations. 

Conclusions. By developing a non-perturbative boso- 
nization approach for bilayer QHSs a,t vt — 1 and by 
considering the fully interacting boson model, we are able 
to properly explore the conjecture of a BEC of excitons 
in the bilayer QHS aX, vt = 1. We show that an ex- 
citon condensate is stable only at very small values of 
the ratio d/l, indicating that the bilayer undergoes a 
phase transition (or crossover) at {d/l)c. A new can- 
didate for the ground state of the system in the region 
{d/l)c < d/l < {d/l)[~^ is requested because its thor- 
ough understanding will be important for the correct de- 
scription of the incompressible-compressible phase tran- 
sition. Details of the above calculations will be presented 
elsewhere. 

We would like to thank H. T. C. Stoof, M. O. Goer- 
big, E. Rezayi, A. Muramatsu, and N. P. Proukakis for 
helpful discussions. R.L.D. and A.O.C. kindly acknowl- 
edge Fundagao de Amparo a Pesquisa do Estado de Sao 
Paulo (Fapesp) and Conselho Nacional de Desenvolvi- 
mento Cientifico e Tecnologico (CNPq) for the financial 
support. 



[1] 

[2] 

[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 
[12] 
[13] 
[14] 
[15] 

[16] 
[17] 
[18] 



For an introduction, see Perspectives in Quantum Hall 
Effects, edited by S. Das Sarma and A. Pinczuk Wiley, 
New York (1997). 

J. P. Eisenstein and A. H. MacDonald, Nature 432, 691 

(2004) ; J. P. Eisenstein, Science 305, 950 (2004). 

S. Q. Murphy et a;.,Phys. Rev. Lett. 72, 728 (1994). 
S. M. Girvin, cond-mat/0108181 

I. B. Spielman et ai, Phys. Rev. Lett. 84, 5808 (2000). 
M. Kellogg et ai, Phys. Rev. Lett. 93, 036801 (2004). 
R. D. Wiersma et ai, Phys. Rev. Lett. 93, 266805 (2004). 
X. G. Wen and A. Zee, Phys. Rev. Lett. 69, 1811 (1992). 
H. A. Fertig, Phys. Rev. B 40, 1087 (1989). 
A. H. MacDonald, P. M. Platzman, and G. S. Boebinger, 
Phys. Rev. Lett. 65, 775 (1990). 

Y. N. Joglekar, and A. H. MacDonald, Phys. Rev. B 64, 
155315 (2001). 

H. A. Fertig and G. Murthy, Phys. Rev. Lett. 95, 156802 

(2005) . 

R. L. Doretto, A. O. Caldeira, and S. M. Girvin, Phys. 
Rev. B 71, 045339 (2005). 

C. Kalhn and B. I. Halperin, Phys. Rev. B 30, 5655 
(1984). 

A. L. Fetter, and J. D. Walecka, Quantum Theory of 
Many-Particle Systems, (Dover Publications, Inc., Mine- 
ola. New York, 2003). 

H. T. C. Stoof and M. Bijlsma, Phys. Rev. E 47, 939 
(1993). 

For a comprehensive review see H. Shi and A. Griffin, 
Phys. Rep. 304, 1 (1998). 

X. M. Chen and J. J. Quinn, Phys. Rev. B 45, 11054 
(1992). 



